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ABSTRACT 

We derive an analytic phenomenological expression that predicts the final mass of the black-hole remnant 
resulting from the merger of a generic binary system of black holes on quasi-circular orbits. Besides recovering 
the correct test-particle limit for extreme mass-ratio binaries, our formula reproduces well the results of all the 
numerical-relativity simulations published so far, both when applied at separations of a few gravitational radii, 
and when applied at separations of tens of thousands of gravitational radii. These validations make our formula 
a useful tool in a variety of contexts ranging from gravitational-wave physics to cosmology. As representative 
examples, we first illustrate how it can be used to decrease the phase error of the effective-one-body waveforms 
during the ringdown phase. Second, we show that, when combined with the recently computed self-force 
correction to the binding energy of nonspinning black-hole binaries, it provides an estimate of the energy 
emitted during the merger and ringdown. Finally, we use it to calculate the energy radiated in gravitational 
waves by massive black-hole binaries as a function of redshift, using different models for the seeds of the 
black-hole population. 

Subject headings: black-hole physics — relativity — gravitational waves — galaxies 



1. INTRODUCTION 

Black-hole (BH) mergers play a central role in 
gravitational-wave (GW) astrophysics, because they are 
expected to be among the main sources for existing and 
future detectors. More specifically, the LIG O/Virgo detec- 
tors ( |Abbott et aT1|2009[ |Acernese et al.|2008) are expected 
to detect mergers of stellar-mass BHs happening within 
several hundred Mpc, when operating in their advanced 
configurati ons. Similarly, future sp ace-based detectors such 
as LISA flAmaro-Seoane et al.|2012| l or DECIGO ( |Kawamura| 
|et al.|2011) will detect mergers of massive BHs (MBHs) up 
to redshifts as high as z ~ 10 or beyond. Even intermediate- 
mass BHs (IMBHs), provided they exist, will be within reach 
of GW detectors, e.g. IMBH-MBH binaries will be detectable 
by LISA or DECIGO, while LMBH-IMBH binaries will be 
detectable with DECIGO or with the planned ground-based 



Einstein Telescope (Punturo et al. 2012, Sathyaprakash et al. 
|20l2) . 

Given their relevance for GW astrophysics, it is not sur- 
prising that BH binaries have received widespread attention 
over the past few years. Because a detailed understanding 
of the dynamics of these systems is crucial in order to pre- 
dict accurately the gravitational waveforms, which, in turn, 
is necessary to detect the signal and extract information on 
the physical parameters of the binaries, numerical simulations 
have been performed by a number of groups for a v ariety of 
mass ra tios, BH spin magnitudes and orientations [see Pfeiffer 
[|( |2012| l for a recent review]. 

However, even today, numerical-relativity simulations are 
computationally very expensive and not able to cover the full 
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seven-dimensional space of parameters of quasi-circular BH 
binaries. Fortunately, phenomenological models have been 
very successful at reproducing many aspects of the dynamics 
of BH binaries as revealed by the numerical simu lations. For 
instance, hybrid "pheno menol ogical waveforms" (Ajith et al. 
2008; Santamaria et al. 2010), i.e., templates that represent 
phenomenological combinations of Post-Newtonian (PN) and 
numerical-relativity (NR) waveforms, can reproduce with 
high precision the NR waveforms for a wide range of binary 
parameters. Similar results are achieved by the effective-one- 
body (EOB) model, which attempts to reproduce not only the 
gravitational waveforms, but also the full dynamics of BH bi- 
naries during the inspiral, mer ger and ringdown phases, by 



resumming the PN dynamics ( Buonanno & Damour| 1999 
Damour et al. 2009 ), and more recently the self-force dynam- 
ics (Baraus se et al.|2012) . 

Other aspects of the dynamics of BH binaries have been 
phenomenologically understood by using combinations of 
PN theory, symmetry arguments, as well as hints from the 
test-particle limit and fits to numerical simulations. For in- 
stance, the final spin magnitude of the BH remnant can be 
predicted by a number of phenomenological formulas (|Rez- 
|zolla et all2008a|b[ |Tichy & Marronetti|[2008l |Buonanno et 
|al.||2008| |Kesden[|2008[ |Rezzolla et al .||2008cl |Barausse & 



Rezzolla 2009 ), starting from the configuration of the binary 
either at small separations r < 10M, or at large separa- 
tions^]?' ~ 10 4 M. These formulas also predict the orientation 



2 For MBHs, the latter are roughly the separations at which the dynamics 
starts being dominated by GW emission, and represent therefore the separa- 
tions at which these phenomenological formulas should work in order to be 
useful in cosmological contexts. 
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of the final spin with good accuracy when applied to small- 



separation binaries, while the formula of Barausse & Rezzolla 
(2009) is also accurate when the binary has a large separa tion, 
e.g. r ~ 10 4 M, in a large portion of the param eter space (Ba- 
rausse & Rezzolla||2009| |Kesden et al.||2010| >. Similar pEe 1 
nomenological formulas have also been proposed for the re- 
coil imparted to t he final BH remnant from the anisotropic 
emission of GW s ( Herrm annet al.|2007bl|Koppitz et al.|2007t 
Rezzolla et al. 2008b ; Campanelli et al. 2007a; Gonzalez et 
al.|2007||Ca mpanelli et al. 200 7b| |Lousto%"2 lochowe r|2009[ 
201 1| |Baker et al.|2007| |2008a| |van Meter et al.||2010| . Be- 
cause most of the anisotropic GW emission takes place as a 
result of the strongly nonlinear merger dynamics, these recoil 
formulas are not predictive, as they depend on quantities that 
can only be derived with full NR simulations, but they are still 
useful in the stat istical studies usually performed in a cosm o- 



logical context (Barauss e|2012| Lousto et a l.|2010b| |2012). 

The dependence of the final mass of the BH remnant on 
the binary's initial paramete rs has also been investigated sys- 
temati cally in the literatur e (Tichy & Marronetti 2008; Boyle| 
& Kesden|[2 008 ; Rei sswig et al.|2009[ |Kesden|2008[ |Lousto| 
et al.| |2010ap] but the knowledge of thi s de pendence is far 
ess detailed. For instance, the formula of Tichy & Marronetti 
(2008) [who built upon previous work by Boyle & Kesden 
(2008 )] is calibrated to reproduce NR results for comparable- 
mass binaries, but does not have the correct test-particle limit 
and is therefore inaccurate for binaries with small mass ratios. 
The formula of Kesden ( 2008 ), on the contrary, has the correct 
test-particle limit, but does not reproduce accurately the NR 
results for comparabl e-mass binaries. Finally, the formula of 
Lousto et al.|([2010a) depends, for generic binary configura- 



tions, on quantities that can only be calculated using full NR 
simulations, and is therefore only useful in statistical studies. 

We here introduce a new phenomenological formula for the 
final mass of the BH remnant (Section^, which, by construc- 
tion, reproduces both the test-particle limit and the regime of 
binaries with comparable masses and aligned or antialigned 
spins, which has been extensively investigated by NR calcu- 
lations. In Section [3] we show that this novel formula repro- 
duces accurately allof the available NR data (even for generic 
spin orientations and mass ratios), both when applied to small- 
and large-separation binary configurations. Furthermore, in 
Section |4] we consider three different areas where our for- 
mula can be useful: (i) we show that it can help reduce the 
phase error of the EOB waveform s during the ringdown; ( ii) 
we combine it with the results of |Le Tiec et al.| ( |2012[ ) for the 
self-force correction to the binding energy of nonspinning BH 
binaries and derive an estimate for the energy emitted during 
the merger and ringdown by nonspinning binaries; (Hi) us- 
ing a semi-analytical galaxy-formation model to follow the 
coevolution of MBHs and their host galaxies, we use our for- 
mula to predict the energy emitted in GWs by MBH binaries 
as a function of redshift, and show that these predictions are 
strongly dependent on the model for the seeds of the MBH 
population at high redshifts. Our final conclusions are drawn 
in Section [5] 

Throughout this paper, geometrized units G = c = 1 are 
used. 

2. THE DEPENDENCE OF THE FINAL MASS ON THE SPINS AND 
THE MASS RATIO 

3 An initi al expre ssion for the radiated energy was also suggested by |Buo-| 
nanno et al. (2007a I, but was restricted to nonspinning binaries and based on 
early NK calculations. 



When deriving a simple algebraic formula that expresses, 
with a given precision, the mass/energy radiated by a binary 
system of BHs, two regimes are particularly well-understood. 
On the analytic side, in fact, the test-particle limit yields pre- 
dictions that are well-known and simple to derive. On the 
numerical side, the simulations of binaries with equal-masses 
and spins aligned or antialigned with the orbital angular mo- 
mentum are comparatively simpler to study, and have been 
explored extensively over the last few years. Hence, it is nat- 
ural that any attempt to derive an improved expression for the 
radiated energy should try and match both of these regimes. 
This is indeed what our formula will be built to do. 

Let us therefore start by considering the test-particle limit 
and, in particular, a Kerr spacetime with mass mi and spin pa- 
rameter a = Si /m\, and a particle (or small BH) with mass 
77X2 on a equatorial circular orbit with radius r ^> To 
first approximation (i.e., for mass ratios q = 7712/ mi <C 1), 
the particle will inspiral towards the BH under the effect of 
GW emission, moving slowly ("adiabat ically") through a se- 
quence of equatorial circular orbits (Kennefick & Ori 1996 ) 
until it reaches the innermost stable circular orbit (ISCO), 
where it starts plunging, eventually crossing the horizon. The 
energy i? la d emitted by the particle during the inspiral from 
r mi to the moment it merges with the central BH can be 
written as 



A I 



= 1 



(1) 



(2) 



f£ co (a) = 3 + Z 2 - sign(a) y/(3 - Zi)(3 + ~Z X + 2Z 2 ) . 



Z x = l + (l-a 2 ) 1/3 (l + a) 1/a + (l-a) 



a/3 



a/3 



3a 2 
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(3) 



(4) 
(5) 



Here, M = mi + m 2 is the total mass, v = m\m 2 j M 2 is 



the symmetric mass ratio, E ISCO and 



are respectively 



the energy per unit mass at the ISCO and the ISCO radius 
in units of mi (Ba rdeen et al.|[l972) , while the remainder, 
o(v), contains the higher-order corrections to the radiated en- 
er 8>0 These corrections account, for instance, for the conser- 
vative s elf-force effects, which affect th e ISCO position and 
energy ( |Barack & Sago 2009 ; |Le Tiec et aT7||2012) >, but also 
for the deviations from adiabaticity, which arise because of 
the finiteness of the mass m2 and which blur the sharp transi- 



tion between inspiral and plunge fBuonanno & Damour 2000 
|Ori & Thorne|2000HKesden|2011| , and, more in general, for 



the energy emitted during the plunge and merger (Berti et al. 
[20071 |Buonanno et al.|2007a|b) . 

If the particle is initially on an inclined (i.e., non-equatorial) 
circular orbit, GW emission will still cause it to adiab atically 
inspiral th rough a sequence of circular orbits (Kennefick & 
|Ori|1996) . Also, the inclination of these orbits relative to the 

4 Without loss of generality, we can assume that the particle moves on a 
prograde orbit (i.e. in the positive-^ direction), and let the spin of the Kerr 
BH point up (a > 0) or down (a < 0). 

5 We here use the Landau symbol o, so that / = o(g) indicates that 
fig — > when g — > 0. Similarly, we will also use the Landau symbol 
O, where instead / = O(g) indicates that f/g — > const when g — > 0. 
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equatorial plane, which can be defined as (Hughes 2001 



cos(i) = 



(6) 



with Q and L z being respectively the Carter constant and 
the azimuthal angular moment um, will remain approximately 
constant during the inspiral (Hughes 2001; Barausse et al. 
2007| l. As in the equatorial case, the particle plunges when 
it reaches the ISCO corresponding to its inclination i. Un- 
like in the equatorial case, though, the radius of the ISCO as a 
function of a and i can only be found numerically. An analyt- 
ical expression, however, can be derived if one considers only 
the spin-orbit coupling of the particle to the Kerr BH, i.e., if 
one considers small spins a <C 1. In that case, in fact, one can 
explicitly check [u sing, for instance, equations (4)-(5) of |Ba-| 
rausse et al. ( 2007[ )] that the ISCO location and energy depend 
only on the combination acos(t), so that at 0(a) 2 , the gen- 
eralization of expressions ([T]i-(|5]l to inclined orbits is given 
by 



Er. 



ad 



71/ 



E„ 



, (a, l)] v + o(v) , 



^ISCo(«: 



'1- 



3r ISCO (a,t) 



(7) 



(8) 



^iscoK ~ r* q co (acos(i)) : 



(9) 

where is given by Q. Expressions (|7]i-(|9) reduce to 

equations (TU-ipj) in the case of equatorial orbits (l = 0) and 
are therefore exact in that limit, with the exception of the 
higher-order terms in v. 

As mentioned above, another case in which we know ac- 
curately the total energy emitted in GWs is given by binaries 
of BHs with equal masses and spins aligned or anti aligned 
with the orbital angular momentum. Reisswig et al. (2009), 
for instance, showed that the energy emitted by these bina- 
ries during their inspiral (from infinite separation), merger and 
ringdown can be well described by a polynomial fit 



M 



= p +pi(ai + a 2 ) +V2{a\ + a 2 ) 2 , (10) 



where a± and a 2 are the projections of the spin parameters 
along the direction L of the orbital angular momentum (cij 
is therefore respectively positive/negative when the spin is 
aligned/antialigned with L), and where the fitting coefficients 
were found to be (Reissw Fg et al.|20 09) po = 0.04826, p% = 
0.0155 9 and p 2 = 0.00485, w ith uncertainties on the order of 
~ 5% (Re isswig et al.||2009| >. We recall that the coefficient 
po can be interpreted as the nonspinning orbital contribution 
to the energy loss (which is the largest one and ~ 50% of the 
largest possible mass loss, which happens for a\ = a 2 = 1), 
Pi can instead be interpreted as the spin-orbit contribution 
(which is < 30% of the largest possible loss), while p 2 can be 
associated to the spin-spin contribution (which is < 20% of 



the largest possi ble loss). Although the fit proposed by Reis 
swig et aT] ( |2009| l predicts a (shallow) minimum for the ra 
diated energy i/ r ad at (ai + o 2 )/2 ~ —0.8, this minimum 
is (very likely) just a n artifact of the fit d ue to the scarce data 
available at that time (Reisswig et al. 2009} . Having now more 
data to analyze, we can enforce the monotonicity of E ra a as a 

6 As in the equatorial case, we can consider only prograde orbits (0 < t < 
it/2) and allow a to be either positive or negative. 
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FIG. 1. — Top panel: Radiated energy, £ ra< j/Af, as a function of the to- 
tal spin of the system along the orbital angular momentum, \ai \ cos/3 + 
| a,2 1 cos 7, for all published NR simulations with q = 1, both with 
aligned/antialigned spins (in red) and for misaligned spins (in blue ). Shown 
instead with a black solid line is the prediction of expression j!4( with the 
coefficients fitted from aligned/antialigned binaries. Bottom panels: resid- 
uals of the NR data from the fitting expression and the corresponding error 
relative to -E ra d /M. 



function of a\ + a 2 by assuming p 2 = 
the minimum of E ISi d to be at (ai 4 
ingly, a fitting expression of the type 



: Pi /4, which constrains 
a 2 )/2 = —1. Interest- 



E r 



ad 



M 



Po +px(a\ 



0-2) 



Pi 
4 



ai + a 2 ) 



(11) 



provides an estimate of the radiated energy which is as accu- 
rate as the one obtained with (JTOj. Indeed, with fitting param- 
eters 

Po = 0.04827 ±0.00039, Pl = 0.01707± 0.00032 , (12) 

this expression reproduces all of the available NR data^j for 
the energy emitted by equal-mass binaries with aligned or an- 
tialigned spins, to within ~ 0.005A/ (except for almost maxi- 
mal spins, see below). Such an accuracy is comparable to the 
typical accuracy of the data themselves, so we can conclude 
that expressions ([TT|i-( [T2"ji summarize our complete knowl- 
edge of the GW emissionTrom this class of binaries to date. 

We note, however, that higher-order terms in the spins may 
be needed in equation (JTTJ to reproduce the data for nearly 
extremal spins. In fact, the maximum value for the radiated 
energy predicted by our fit, i.e., 9.95% of the total mass of 
the binary at infinite separation when ai = a 2 = 1, is signif- 



icantly less than the 10.95% found by Lovelace et al. (2012 1 
for a\ = a 2 ~ 0.97. Such a large value for i-' ra d is some- 
what off the general trend shown by the other NR data for 
large aligned spins. However, it is clear that higher-order spin 

7 The NR data considered is relativ e to the following references listed in 
alphabetical o rder: |Baker et al.|( 2008b i; Berti et al. 1 2007 2008 1; Campanelli 
letaLl (20061; | Chu et al.N!>0Wt; |lJhu]l|2012t; |Mannam et al.l|20uij|fi01uF 
|Kelly et al.M201 1>; ILovelace et al.M201 1| |2 012 1; Marronett ieTaTT^OOijt ; 
|Pollney et al.|j2W/} ; |Pollney & Reisswig|{204 TJ ; |Reisswig et al.||200^ ~~^ 
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terms may have to be added if more numerical data for high- 
spin configurations becomes available and confirm this result. 

Using therefore the knowledge of the radiated en- 
ergy from the test-particle limit and from the equal-mass 
aligned/antialigned configurations, we derive an expression 
valid for generic binaries. As a first step, let us note that the 
PN binding energy of an equal-mass binary of spinning BHs 
depends o n the spins, at 1 .5 PN order , i.e., at leading order in 
the spins (Bark er"& Q'Conneil|1975| l, only through the com- 
bination 



L-(Si + S 2 ) |tti| cos/3 + \a 2 \ cos 7 
M 2 ~ 4 



(13) 



where |ai| and \a 2 \ are the spin magnitudes, and /3, 7 are the 
angles between the orbital angular momentum unit vector L 
and the spins of the first and second BH, respectively. We can 
therefore attempt to extend expression ( 1 1 1 to generic equal- 
mass binaries simply by replacing a\ + a 2 with |oi| cos /3 + 
\a 2 \ cos 7, i.e., 



-Brad 



= Pa +pi(|ai| cos/3 + \a 2 \ cos 7) 



Pi 



|ai| cos/3 + \a 2 \ cos^Y 



(14) 



As a check of this ansatz, in the top panel of Fig. [T] we have 
plotted the radiated energy, E ia ^/M, as a function of the to- 
tal spin along the orbital angular momentum, |ai|cos/3 + 
|a.2|cos7, for all published NR simulations with q = 1, 
both with aligned/antialigned spins (in red; see footnote [7]i 
and with misaligned spins (in bkuQ- Also, we show with 
a black solid line the prediction of expression ( fl4) with the 
coefficients fitted from aligned/antialigned binaries [equation 



12 1], In the bottom panels we show instead the residuals of 
theNR data from the same curve and the corresponding errors 
relative to _E ra( j/M. Clearly, while future simulations that 
are more accurate or describe more involved configurations 
may present deviations from our simple ansatz, all published 
simulations for equal-mass binaries are in reasonable agree- 
ment with equation ( 14 1, with residuals of < 1% and errors 
of < 10% relative to trie radiated mass. Note that these errors 
are comparable with the intrinsic scatter of the different NR 
data. 

Because in the test-particle limit the angle /3 becomes the 
angle between the spin Si of the Kerr BH and the orbital 
angular momentum of the particle, thus coinciding with the 
angle t defined in (|6j, it is natural to rewrite equations (|7]l-((9]l 



E 



rad 



M 



I 1 - Asco («)] v + °0) 



£isco(a) = a/1 - 



ISCO v ' 

where we have defined 

_ _ L ■ (Si + S 2 ) _ \ai\ cos/3 + q 2 |a 2 | COS7 



M 2 



(i + # 



(15) 



(16) 



(17) 



If we now assume that the higher-order term o(y ) in equation 
( 15 1 is quadratic in v, we can determine it by imposing that 



8 The NR data for equal-mass misaligned bi naries are relative to th e 
following references listed in alphabetical order: [Herrmann et al.| 1 2007a I ; 
|Lousto et : al.]{2^T2) ; [TTc"hy & M arronetti (2007 2 00^ 



we recover the equal-mass expression ( 14 1 for q 
obtaining the final expression 



1, thus 



E r 



ad 



M 



,(a)]u 



+ 4 v 2 [4p + 16p!o(o + 1) + E ISCO (a) - 1] , 

(18) 

where E JSCO (a) is given by (TTSJ). By construction, there- 
fore, expression ([18) has the correct behavior both in the test- 
particle limit and for equal-mass binaries. Also, we stress that 
the fitting coefficients [given by ( [T2"] >] are obtained using only 
a subset of the NR data (i.e., those for equal-mass binaries 
with aligned/antialigned spins). 

3. COMPARISON TO DATA: BINARIES AT SMALL AND LARGE 
SEPARATIONS 



In order to test the accuracy of expression ( 18 I, we used the 
data of 186 numerical simulations of inspiraTling and merg- 
ing BH binaries^] which have reported the ratio Mf/M 
I — E Ta d/M between the final mass of the BH remnant, Mf, 
and the mass M = mi + m 2 of the binary at infinite separa- 
tion. In cases where this ratio was not reported explicitly, we 
have reconstructed it from the energy radiated during the nu- 
merical simulation using PN expressions. More specifically, 
we calculate the radiated energy as 



#rad = E Iad + l^bind(^o)| , 



(19) 



where E hind (fl) is the 3PN binding energy as function of 



the orbital frequency 51 (Buonanno et al. 2003b), and fin. is 
the initial orbital frequency of the simulation (either reported 
explicitly or, when unavailable, reconstructed from the ini- 
tial puncture data). In addition, in those cases where the 
simulation results are normalized in terms of the Arnowitt- 
Deser-Misner mass Majjm, we approximate it as Madm = 



For each binary, we apply our expression ( 18 1 to the ini- 
tial configuration of the numerical simulation (where the bi- 
nary typically has a " small separation" r < 10M). How- 
ever, in order to check whether our expression predicts the 
final mass correctly also for widely separated binaries, we 
have also integrated the initial binary back to a "large sepa- 
ration" r = 2 x 10 4 M using the quasi-circular PN evolu- 
tion equations of Buonan no et al.| ( |2003b| l (which are accu- 
rate through 3.5PN order for the adiabatic evolution of the 
orbital frequency, and through 2PN order for the dynamics of 
the spins). For massive BHs, this is roughly the separation at 
which the dynamics starts being dominated by GW emission, 
and is therefore the separation at which our expression ( 18 1 
ought to w ork if we want it to be useful in cosmological con- 
texts [cf. Barausse & Rezzolla (2009) and the discussion in 
|Barausse|p312b ]. 

The results of these comparisons are summarized in Pig. [2j 
which shows the difference between the data and our expres- 
sion, both for small (left panel) and large separations (right 

9 The data is relative t o the following references listed in alphabetical or - 
der|Baker et al.|j20 08bl: Berti et al. (2008 2007 1; Buchman eTaTl <20T2t; 
Campan eHTet al.H2 006 2009 1; ; Chu et al.M20U9»;|ChuH2012»;|Oonzalez e t 
al l jiMjty ; IMannam et al.j |2W8l I^OlOt; |Herrmann et al N20U7at; |kelly et 
[ aT*1(2UI l(;|Lousto & Zlochower|l2UU9);|Lousto et al.H2U12);|Lovelace et al. 
)11| 12012); IMarronetti et al.|l| 2008 l; Na kano et al.M20lUnFollney et al. 
nn7)nPoilney & ReisswigH2011) ; |keisswig et al.||2U09( ; [Tichy & Mar- 
|ronem|f20U7||l2U0yi. 
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FIG. 2. — Residuals for the fitting formula at small and large separations as a function of a dummy index representing the binaries in our dataset. Binaries 
with spins aligned/antialigned with the orbital angular momentum are plotted in red, while binaries with misaligned spins are shown in blue. Cyan and violet 
lines represent the la and 2<r errors (estimated a posteriori) of the data with spins aligned/antialigned with the orbital angular momentum. The first 50 points 
correspond to the simulations performed after 2010. 



panel) as a function of a dummy index representing the or- 
dering of the binaries in our dataset. As can be seen the re- 
sults at small and large separations are almost indistinguish- 
able. This does not come as a surprise, because the projection 
of the total spin on the direction of the angular momentum 
[equation ( fT3| l] is approximately conserved during the inspi- 
ral, in most of the parameter space of quasi-circular binaries 
(see discussion in Barausse & Rezzolla (2009 ) for more de- 
tails). Binaries with spins aligned/antialigned with the orbital 
angular momentum are plotted in red, while binaries with mis- 
aligned spins are shown in blue. Also shown are cyan and vi- 
olet lines representing the la and 2a errors of the data with 
spins aligned/antialigned with the orbital angular momentum, 
as obtained a posteriori by comparing them to the fit (jTTJ (cf. 
Fig. [TJ. Also, the first 50 points correspond to the simula- 
tions performed after 2010, while others correspond to the 
simulations performed in 2006-2009. Although the quality 
of numerical simulations improved substantially in last few 
years, the "old" data give residuals comparable to the "new" 
ones. Furthermore, the residuals for the binaries with spins 
aligned/antialigned with the orbital angular momentum ap- 
pear to be comparable with those for the binaries with mis- 
aligned spins. 

We stress that while our expression is in reasonable agree- 
ment with all the published data, both at large and small sep- 
arations, there are still large gaps in the parameter space of 
BH binaries that prevent us from testing our approach more 
thoroughly. This is best seen in Fig. [3] where we plot the fi- 
nal mass of the remnant for all the published data for binaries 
with a\ cos (3 — a-i cos 7 (blue circles), as well as the predic- 
tions of our expression when applied to the "small-separation" 
initial data of the simulations (meshed surface). Clearly, spin- 
ning binaries with unequal mass ratios are essentially absent, 
and simulations for such binaries will provide a very signifi- 
cant check of our expression ( 18 1. Nevertheless, the simple 
functional dependence shown by the available data, whose 



behaviour can be well captured with low-order polynomials 
(with the possible exception, as we already stressed, of almost 
maximally spinning configurations), is quite remarkable. 
The graphical representation of the data in Fig. [3] also al 



lows t o reinforce a remark already made by |Reisswig et al. 
( 2009 1, namely, that the largest radiated energy, i/' r ad(o = 
'. ~j]M — 9.95%, is lost by binaries with equal-mass and max- 
imally spinning BHs with spins aligned with the orbital an- 
gular momentum. Hence, BH binaries on quasi-circular or- 
bits are among the most efficient sources of energy in the 
Universe. Note, however, that equal-mass binaries are not 
always the systems that lose the largest amount of energy. 
Indeed, unequal-mass systems with sufficiently large spins 
aligned with the angular momentum can lead to emissions 
larger than those from equal-mass binaries but with large an- 
tialigned spins. For instance, a binary with v = 0.15 and 
di = a,2 = 1 will radiate more than a binary with v = 0.25 
and a\ = —ai- 

Notwithstanding the limited coverage of the parameter 
space, we can note that our approach substantially im proves 



upon earlier fo rmulas for the final mass. For instance, Tichy 
& Marronetti ( 2008 ) [building on the work of Boyle & Kes 



|den| (|2008 )] suggested a formula linear in the symmetric mass 
ratio v, but the coefficients needed to fit NR results are such 
that the test-particle limit ([lJ-Q is not recovered. As men- 
tioned earlier, because we recover the test-particle limit ex- 
actly, our expression reproduces the published data more ac- 
curately, especially for small mass-ratio configurations (cf. 
the discussion on the effective-one-body model in t he next 
section). An other example is given by the formula of |Lousto] 
|et al.| ( |2010a[ ), which has the correct test-particle limit but de- 
pends, for generic configurations, on parameters that describe 
the binary's plunge and merger and which can only be cal- 
culated with numerical simulations. Our algebraic formula, 
instead, allows one to calculate the final mass with reasonable 
accuracy, using only information on the initial binary config- 



6 



Barausse, Morozova and Rezzolla 




a 1 cosp = a 2 cosy 



FIG. 3.— Mass of the final BH, M 



f 



M 



-E ra d, and con'esponding 
fit for all the published binaries with a± cos (5 = ai cos 7. Note the simple 
functional dependence of the -E ra d> whose behaviour can be well captured 
with low-order polynomials. 



uration, at any separation. 

4. APPLICATIONS OF THE NEW FORMULA 

In the following Sections we discuss three different exam- 
ples of how our new expression for the energy radiated during 
the inspiral, merger and ringdown of two BHs can be used in 
contexts that range from GW physics to cosmology. 



4.1. Merger-ringdown energy 



We can combine our expression (|T8]i f or the total radiated 
energy with the recent results of LeTiec et al. ( 2012[ ) for the 
binding energy of a binary of nonspinning BHs at next-to- 
leading order in the mass ratio, and obtain an expression for 
the energy emitted in the merger and ringdown phases of non- 
spinning BH binarie s. More specifically, for these binaries 
Le Tie c et aL] p012) found the total energy (i.e., the binary's 
mass M at infinite separation, plus the binding energy) to be 



E{x) = M 



( l-2x 



1 )u + u z E SF {x) 



0{uf 
(20) 

where x = (M Q) 2 ^ 3 and fl is the orbital frequency. The 
self-force contribution reads 



E SF (x) 



SF 



(x) 



-l+Vl - 3a: + 



2 3 
where zsf is given by the fitting function 

2x (1 + ciix + CL2X 2 ) 



1 + 032; + a^x 2 + 05X 3 ' 



x{7 - 2Ax) 

6(1 - 3x) 3 / 2 ' 
(21) 



(22) 



which is accurate to within 10 5 with a\ = —2.18522, 
a 2 = 1.05185, a 3 = -2.43395, C14 = 0.400665, and 
a 5 = —5.9991, and where we use a prime to denote deriva- 
tives with respect to x. The minimum of E ( x) marks the lo- 
cation o f the ISCO [see |Le Tiec et al.|(2012| ); |B"uonanno et al.| 
P003a| l] and lies at 



l + ^Cn) +0(v 2 ), (23) 



with Co = 1.2510(2). Replacing x lsco in equation (20 1, one 
obtains that the energy emitted during the inspiral is 



-Erad 



insp 



-^(•^isco ) 




0.037763 v 2 +(D(vf 



0.057191 v + 0.037763 v 2 + 0{v) 



(24) 



Expressing now the total radiated energy as the sum of the 
one lost during the inspiral and the one lost during the plunge- 
merger and ringdown, i.e., 



-^rad — -^rad,insp ~T" -^rad, merger— rd • 



(25) 



and using equation ( 18 1, we obtain an expression for the en- 
ergy radiated during the plunge-merger and ringdown of non- 
spinning BH binaries: 



-Erad .1 



M 



— 0.506 v 2 , 



(26) 



For equal-mass binaries, this energy is almost twice the en- 
ergy lost during the whole inspiral. In expression (j26j we have 
neglected terms of order 0(v) 3 , so in principle this equation 
may not hold for comparable-mass bi naries. However, the 
binding energy ( 20 1 has been found by |Le Tiec et al.| ( |2012| l 
to be in very goodagreement with NR results for comparable 
masses, and we ther efore expect the same to hold for expres- 
sion (j26j>. Indeed, [Berti et al. (2007) have estimated that the 
energy emitted by equal- and unequal-mass nonspinning BH 
binaries after the 3PN ISCO is given by 



-^rad.merger-rd ^ „ y 2 



M 



(27) 



Even more strikingly IBerti et al.| ( |2007| l; |Nagar et al.| (2007); 
Bernuzzi & Nagar| ( |20l0| ) showed that the energy emitted by 
a particle in a Schwarzschild spacetime during its plunge- 
merger is given by 



plunge 

M 



0A7v z 



(28) 



in reasonable agreement with (27i and (26 1. As a result, at 
least in the nonspinning case, one can reproduce the radiated 
energy predicted by our final expression ( 18 1 for comparable- 
mass binaries, simply by summing the energy emitted by a 
particle during the inspiral [equation ( [24] )] and during the 
plunge-merger [equation (|28)]. 

This confirms previous work showing that perturbative re- 
sults, when expressed in terms of the symmetric mass ratio 
v, are in good a greement with NR results, even for compara- 
ble ma s ses [see[Detweiler & Szedenits| ( |1979| l; 



1995 



Anninos et al. 



Berti et al.| ( |2007| l for the GW fluxes, |Le Tiec et al 
for the periastron precession, and |Le Tiec et al.| ( |2012| l 



( 2011 

1 or the binding energy]. 



4.2. 



Tuning of the effective-one-body model 

We recall that the EOB ( |Buonanno & Damour||1999| > is a 
phenomenological model aiming at describing the dynamics 
and waveforms of BH binaries, during the inspiral, merger 
and ringdown phases, combining information from PN theory, 
perturbative calculations and NR. While developed initially 
for nonspinning BH s (Buonanno & Damour 1999 Damour 
et al. 2000; Barausse et al. 2012), the model has been more 

I II r I II 1 



On the mass radiated by coalescing black-hole binaries 



7 



recently generalized to spinning ones (Damour 2001 ; Damour 



et al.|2b08||Barausse & Buonanno|20iO|j2011MNagar]2011| i. 

To accurately describe the ringdown phase, the EOB needs 
expressions predicting the final mass and spin. For non- 
spinning BHs, these could be estim ated self-consistently 
within the EOB along the lines of Damour & Nagar (2007). 
However, a generalization of this approach to spinning BHs, 
especially if precessing, is not straightforward. Moreover, be- 
cause the final BH's mass and spin are used to calculate the 
frequencies and decay times of the quasi-normal modes, even 
small inaccuracies in the prediction of the remnant's mass and 
spin introduce considerable phase errors i n the EOB wave- 
forms. For instance, the EOB model of Tara cchini et al.| 
( |2012 1 was compared to NR waveforms for nonsp inning BHs 
wit h mass rati o q = 1/6, using the formula of Tichy & Mar- 



|ronetti| ( |2008| l to calculate the final mass. Because that for- 
mula does not have the correct test-particle limit [e.g. for non- 
spinning BHs, it predicts E YaA /M = 0.194 v + 0(v) 2 instead 
of E lad /M = 0.057^ + 0{v) 2 ], the EOB waveforms were 
accumulating a phase difference of ~ 0.2 rad from the NR 
ones during the ringdown. With our new formula, this phase 
difference during the ringdown decreases impressively to less 
than 0.05 rad (Taracchini 2012, private communication). 

4.3. GW emission by MBHs 

According to our present understanding of galaxy forma- 
tion, most galactic nuclei should host a massive BH, with 
mass up to 10 10 M©. Information on the masses and spins 
of these BHs can be extracted from p resent electromagnetic 
observations (see e.g. Li et al.^2012) for recent constraints), 
but much more accurate data will be provided by future space- 
based GW detectors such as LISA or DECIGO, or terrestrial 
ones such as the Einstein Telescope. These detectors will 
be able to observe the mergers of MBHs, which take place 
when their host galaxies coalesce. S emi-analytical g alaxy- 
formation models, such as e.g. that of |Barausse] f2012), have 
been employed to understand the MBH merger history and 
therefore the event rates for these detectors, suggesting hun- 
dreds of events per year for DECIGO, from a few to hundreds 
of events per year for LISA, and up to a few tens of events per 
year for the Einstein Telescope. A detector-independent diag- 
nostics of the importance of massive BH mergers, however, is 
given by the energy radiated in GWs by massive BH binaries, 
per unit comoving volume and as a function of cosmic time. 
We have calculated this q uantity using the galaxy-formation 
model of Barausse (2012), our new expression ( flH) , and two 
competing models for the seeds of the massive BHs at high 
redshift - namely a "light-seed" scenario in which the seeds 
have mass M see d ~ 100 M at z = 15 — 20, and a "heavy- 
seed" model in whic h the seeds form a z ~ 15 with mass 
Mseed ~ 10 5 Mq (see Barausse (2012 1 and references therein 
for more details). 

The results for this quantity are shown in Fig. [4] Clearly, 
the heavy-seed scenario predicts much stronger GW emission 
at redshifts z > 3, which is not surprising since mergers hap- 
pen initially between BHs with masses ~ M see d> an d the ra- 
diated energy scales with the total mass of the BH binaries. 
At z ~ 0, instead, the two models yield very similar results, 
because bo th reproduce th e observed local massive BH mass 
function (Barausse 2012). Given the significantly different 
yields in the GW emission expected from these two scenar- 
ios of galaxy formation, future space-borne and terrestrial in- 
terferometers would provide important and unambiguous in- 
formation on the BH population at early redshifts [see also 



o 
ft 

>. 



> 

"3 



S 
T3 



H 

CM 



10 4 



1000 



100 



10 



0.1 



heavy seeds 
light seeds 



10 

z 



15 



20 



FIG. 4. — The energy emitted by massive BH mergers per unit redshift and 
unit comoving volume, as a function of redshift. The two lines refer either 
to the "light-seed" scenario (red solid curve) or to the "heavy-seed" scenario 
(blue dashed line). 



ISesana et al.|p009l|20TT) ; |Gair et al.| ( |2009||201 1| >]. 

Finally, we note that by integrating the results of Fig. HI we 
find that the total energy density in GWs from massive BH bi- 

7.4 x 1O 2 M /Mpc 3 in the 



1.8 x 10 4 M m /Mpc 3 



nailes at z = is p GW , mcrgcrs 

light-seed scenario and pew, mergers ~ ±.o ^ ±u jh 
in the heavy-seed scenario, corresponding to a cosmologi 
cal density parameter f2 G w,raer g ers = PGW,mer ge rs/Pcrit ~ 
5.4 x 10~ 9 (light-seed scenario) or I^gw, mergers ~ 1.3xl0~ 7 
(heavy-seed scenario). 

5. CONCLUSIONS 

We have presented a novel algebraic formula to measure 
the energy radiated by coalescing binary systems of BHs with 
generic spin magnitudes and orientations and arbitrary mass 
ratios. Our expression uses information on the binary con- 
figuration at an arbitrary separation and reproduces correctly 
the two regimes in which the radiated energy is known best, 
namely, the test-particle limit (which is known analytically) 
and the comparable-mass case (which has been extensively 
investigated with NR simulations over the last few years). Be- 
cause it smoothly interpolates these two regimes, we expect 
our formula to work reasonably well also for intermediate 
mass ratios. Indeed, we have verified that it reproduces the 
results of all the NR simulations published so far in the litera- 
ture, including those with unequal masses, to within an error 
which is comparable to the typical errors of the simulations. 
In addition, we have checked that our formula works equally 
well when applied to binaries starting at small separations 
(i.e., r < 10M) and at large separations (i.e., r ~ 10 4 M), 
thus opening up the possibility of using our expression also in 
cosmological contexts. 

The algebraic nature of our expression makes it a useful 
tool in a variety of contexts that range from GW physics to 
cosmology. As representative examples we have discussed 
three different applications, namel y: (i) we hav e shown that, 
when combined with the results of |Le Tiec et al.| ( |2~012| l for the 



8 



Barausse, Morozova and Rezzolla 



self-force correction to the binding energy of nonspinning BH 
binaries, the new formula provides an estimate for the energy 
emitted during the merger and ringdown, and that this esti- 
mate confirms the conjecture that the results of perturbative 
calculations may be successfully extrapolated to comparable- 
mass binaries when expressed in terms of the symmetric mass 
ratio v\ (ii) we have shown that the new formula can help re- 
duce the phase error of the EOB waveforms during the ring- 
down; ( Hi) using a semi-analytical galaxy-formation model to 
follow the coevolution of MBHs and their host galaxies, we 
have used our formula to predict the energy emitted in GWs 
by MBH binaries as a function of redshift, and found that 
these predictions strongly depend on the scenario adopted for 
the MBH seeds at high redshifts, thus making GW emission a 
powerful cosmological diagnostic. 

Additional uses of the new formula can be easily consid- 
ered and a particularly relevant one is the impact of the mass 
loss on the accretion disk surrounding the MBH binary. The 
dynamics of the disk, in fact, can change considerably as a 
result of the very rapid change in the gravitational mass of the 
system, with the formation of large shocks, which ar e poten- 
tially detectable via their electromagnetic emiss ion ( |0'Neill| 
|et al.|2009[ [Rossi et al.|20lU||Zanotti et al.|2010) . 



As a final remark we note that because our approach ex- 
ploits knowledge derived from NR simulations, the accuracy 
of the final-mass formula can be improved as additional and 
more precise NR simulations, especially with highly-spinning 
BHs, become available. 
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viding the final mass o f an unpublished simulation used in 
Taracc hTni et al.| ( |2012| ) and presented in his PhD thesis ( |Chu| 
|2012p . Special thanks go to E. Berti for several useful dis- 
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our formula significantly improves the EOB phase error for 
small mass ratios. E.B. acknowledges support from a CITA 
National Fellowship at the University of Guelph. This work 
was supported in part by the DFG grant SFB/Transregio 7. 



REFERENCES 



Abbott, B. P., Abbott, R., Adhikari, R., et al. 2009, Reports on Progress in 

Physics, 72, 076901 
Acernese, E, Alshourbagy, M., Amico, P., et al. 2008, Classical and 

Quantum Gravity, 25, 184001 
Ajith, P., Babak, S., Chen, Y. et al. 2008, Phys. Rev. D, 77, 104017 
Amaro-Seoane, P., Aoudia, S., Babak, S., et al. 2012, |arXiv:r2 01.3621 
Anninos, P., Price, R. H., Purlin, J., Seidel, E., & Suen, W.-M. 1995, 

Phys. Rev. D, 52, 4462 
Baker, J. G., Boggs, W. D., Centrella, J., et al. 2007, ApJ, 668, 1140 
Baker, J. G, Boggs, W. D., Centrella, J., et al. 2008, ApJ, 682, L29 
Baker, J. G, Boggs, W. D., Centrella, J., Kelly, B. J., McWilliams, S. T., & 

van Meter, J. R. 2008, Phys. Rev. D, 78, 044046 
Barack, L., & Sago, N. 2009, Physical Review Letters, 102, 191101 
Barausse, E. 2012, MNRAS, 423, 2533 

Barausse, E., Buonanno, A., & Le Tiec, A. 2012, Phys. Rev. D, 85, 064010 

Barausse, E., Hughes, S. A., & Rezzolla, L. 2007, Phys. Rev. D, 76, 044007 

Barausse, E., & Buonanno, A. 2010, Phys. Rev. D, 81, 084024 

Barausse, E., & Buonanno, A. 2011, Phys. Rev. D, 84, 104027 

Barausse, E., & Rezzolla, L. 2009, ApJ, 704, L40 

Bardeen, J. M., Press, W. H., & Teukolsky, S. A. 1972, ApJ, 178, 347 

Barker, B. M., & O'Connell, R. F. 1975, Phys. Rev. D, 12, 329 

Bernuzzi, S., & Nagar, A. 2010, Phys. Rev. D, 81, 084056 

Berti, E., Cardoso, V., Gonzalez, J. A., et al. 2007, Phys. Rev. D, 76, 064034 

Berti, E., Cardoso, V., Gonzalez, J. A., Sperhake, U., & Brugmann, B. 2008, 

Class. Quantum Grav., 25, 114035 
Boyle, L., & Kesden, M. 2008, Phys. Rev. D, 78, 024017 
Buchm an, L. T„ Pfe iffer, H. P., Scheel, M. A., & Szilagyi, B. 2012, 

arXiv:1206.3015>'l 
Buonanno, A., Pan, Y., Baker, J. G, et al. 2007, Phys. Rev. D, 76, 104049 
Buonanno, A., Kidder, L. E., & Lehner, L. 2008, Phys. Rev. D, 77, 026004 
Buonanno, A., & Damour, T. 1999, Phys. Rev. D, 59, 084006 
Buonanno, A., & Damour, T. 2000, Phys. Rev. D, 62, 064015 
Buonanno, A., Chen, Y, & Vallisneri, M. 2003, Phys. Rev. D, 67, 024016; 

Erratum: 2006, Phys. Rev. D, 74, 029903(E) 
Buonanno, A., Chen, Y, & Vallisneri, M. 2003, Phys. Rev. D, 67, 104025; 

Erratum 2006, Phys. Rev. D, 74, 029904 
Buonanno, A., Cook, G. B., & Pretorius, F. 2007, Phys. Rev. D, 75, 124018 
Campanelli, M., Lousto, C. O., & Zlochower, Y. 2006, Phys. Rev. D, 74, 

041501 

Campanelli, M., Lousto, C, Zlochower, Y, & Merritt, D. 2007, ApJ, 659, L5 
Campanelli, M., Lousto, C. O., Zlochower, Y, & Merritt, D. 2007, 

Phys. Rev. Lett., 98, 231102 
Campanelli, M., Lousto, C. O., Nakano, H., & Zlochower, Y. 2009, 

Phys. Rev. D, 79, 084010 
Chu, T., Pfeiffer, H. P., & Scheel, M. A. 2009, Phys. Rev. D, 80, 124051 



Chu, T. 2012, Numerical Simulations of Black-Hole Spacetimes, PhD 

thesis, California Institute of Technology 
Damour, T., Jaranowski, P., & Schafer, G. 2000, Phys. Rev. D, 62, 08401 1 
Damour, T. 2001, Phys. Rev. D, 64, 124013 
Damour, T., & Nagar, A. 2007, Phys. Rev. D, 76, 044003 
Damour, T., Jaranowski, P., & Schafer, G. 2008, Phys. Rev. D, 78, 024009 
Damour, T., Iyer, B. R., & Nagar, A. 2009, Phys. Rev. D, 79, 064004 
Detweiler, S. L., & Szedenits, E., Jr. 1979, ApJ, 231, 211 
Gair, J. R., Mandel, I., Sesana, A., & Vecchio, A. 2009, Classical and 

Quantum Gravity, 26, 204009 
Gair, J. R., Sesana, A., Berti, E., & Volonteri, M. 201 1, Classical and 

Quantum Gravity, 28, 094018 
Gonzalez, J. A., Brugmann, B., Sperhake, U., Hannam, M., Husa, S. 2007, 

Phys. Rev. Lett., 98, 231101 
Gonzalez, J. A., Sperhake, U., & Brugmann, B. 2009, Phys. Rev. D, 79, 

124006 

Hannam, M., Husa, S., Brugmann, B., & Gopakumar, A. 2008, 

Phys. Rev. D, 78, 104007 
Hannam, M., Husa, S., Ohme, F, Miiller, D., & Brugmann, B. 2010, 

Phys. Rev. D, 82, 124008 
Herrmann, F, Hinder, I., Shoemaker, D. M., Laguna, P., & Matzner, R. A. 

2007, Phys. Rev. D, 76, 084032 
Herrmann, F, Hinder, I., Shoemaker, D. M., Laguna, P., & Matzner, R. A., 

2007, ApJ, 661, 430 
Hughes, S. A. 2001, Phys. Rev. D, 64, 064004 

Kawamura, S., Ando, M., Seto, N, et al. 2011, Classical and Quantum 

Gravity, 28, 094011 
Kelly, B. J., Baker, J. G, Boggs, W. D., McWilliams, S. T., & Centrella, J. 

201 1, Phys. Rev. D, 84, 084009 
Kennefick, D., & Ori, A. 1996, Phys. Rev. D, 53, 4319 
Kesden, M. 2008, Phys. Rev. D, 78, 084030 

Kesden, M., Sperhake, U., & Berti, E. 2010, Phys. Rev. D, 81, 084054 
Kesden, M. 201 1, Phys. Rev. D, 83, 10401 1 

Koppitz, M., Pollney, D., Reisswig, C. et al. 2007, Phys. Rev. Lett., 99, 
041102 

Le Tiec, A., Barausse, E., & Buonanno, A. 2012, Physical Review Letters, 
108,131103 

Le Tiec, A., Mroue, A. H., Barack, L., et al. 201 1, Physical Review Letters, 
107,141101 

Li, Y.-R., Wang, J.-M., & Ho, L. C. 2012, ApJ, 749, 187 

Lousto, C. O., & Zlochower, Y. 2009, Phys. Rev. D, 79, 064018 

Lousto, C. O., Campanelli, M., Zlochower, Y, & Nakano, H. 2010, Class. 

Quantum Grav., 27, 114006 
Lousto, C. O., Nakano, H., Zlochower, Y, & Campanelli, M. 2010, 

Phys. Rev. D, 82, 129902 (E) 
Lousto, C. O., & Zlochower, Y. 2011, Phys. Rev. Lett., 107, 231102 



On the mass radiated by coalescing black-hole binaries 



9 



Lousto, C. O., Zlochower, Y., Dotti, M, & Volonteri, M. 2012, Phys. Rev. D, 
85, 084015 

Lovelace, G., Scheel, M. A., & Szilagyi, B. 2011, Phys. Rev. D, 83, 024010 
Lovelace, G., Boyle, M., Scheel, M. A., & Szilagyi, B. 2012, Class. 

Quantum Grav., 29, 045003 
Marronetti, P., Tichy, W., Briigmann, B., Gonzalez, J., & Sperhake, U. 2008, 

Phys. Rev. D, 77, 064010 
Nagar, A., Damour, T., & Tartaglia, A. 2007, Classical and Quantum 

Gravity, 24, 109 
Nagar, A. 201 1, Phys. Rev. D, 84, 084028 
Nakano, H., Zlochower, Y., Lousto, C. O., Campanelli, M. 201 1, 

Phys. Rev. D, 84, 124006 
O'Neill, S. M., Miller, M. C, Bogdanovic, T., Reynolds, C. S., & 

Schnittman, J. D. 2009, Astrophys. J., 700, 859 
Ori, A., & Thorne, K. S. 2000, Phys. Rev. D, 62, 124022 
Pollney, D., Reisswig, C, Rezzolla, L., et al. 2007, Phys. Rev. D, 76, 124002 
Pollney, D., & Reisswig, C. 201 1, ApJ, 732, L13 
Pfeiffer, H. 2012, Class. Quantum Grav., 29, 124004 
Punturo, M B., Abernathy, M., Acernese, F., et al. 2010, Class. Quantum 

Grav., 27, 084007 



Reisswig, C, Husa, S., Rezzolla, L., et al. 2009, Phys. Rev. D, 80, 124026 
Rezzolla, L., Diener, P., Dorband, E. N, et al. 2008, ApJ, 674, L29 
Rezzolla, L., Dorband, E. N., Reisswig, C, et al. 2008, ApJ, 679, 1422 
Rezzolla, L., Barausse, E., Dorband, E. N., et al. 2008, Phys. Rev. D, 78, 
044002 

Rossi, E. M., Lodato, G, Armitage, P. J., Pringle, J. E., & King, A. R. 2010, 

Mon. Not. R. Astron. Soc, 401, 2021 
Santamarfa, L., Ohme, E, Ajith, P., et al. 2010, Phys. Rev. D, 82, 064016 
Sathyaprakash, B., Abernathy, M., Acernese, E, et al. 2012, Class. Quantum 

Grav. 29, 124013 
Sesana, A., Gair, J., Mandel, I., & Vecchio, A. 2009, ApJ, 698, L129 
Sesana, A., Gair, J., Berti, E., & Volonteri, M. 201 1, Phys. Rev. D, 83, 

044036 

Tichy, W., & Marronetti, P. 2007, Phys. Rev. D, 76, 061502 

Tichy, W., & Marronetti, P. 2008, Phys. Rev. D, 78, 081501 

Taracchini, A., Pan, Y, Buonanno, A., et al. 2012, Phys. Rev. D, 86, 02401 1 

van Meter, J. R., Miller, M. C, Baker, J. G., Boggs, W. D., & Kelly, B. J. 

2010, ApJ, 719, 1427 
Zanotti, O., Rezzolla, L., & Del Zanna, & Palenzuela, C, 2010, Astron. 

Astrophys., 523, A8 



